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Abstract 



The hierarchy of long-distance matrix elements (MEs) for quarkonium pro- 
duction depends on their scaling with the velocity v of the heavy quark in 

F^ ■ the bound state. Ranges for the velocities in various bound states and un- 

certainties of colour- singlet MEs are estimated in a quark-potential model. 

KA , Different possibilities for the scaling with v of the MEs are discussed; they 

H \ depend on the actual values of v and the QCD scale. As an apphcation, J/^ 

polarization in e"'"e~ annihilation is discussed. The first non-perturbative 
estimates of colour-octet MEs are presented and compared with phenomeno- 
logical determinations. Finally, various predictions of prompt quarkonium 
production at LEP are compared. 
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1 Introduction 

Despite the appearance of genuine non-perturbative elements, we can calculate the pro- 
duction of a heavy quarkonium bound state H{J^^): the cross section is given as a sum 
of terms 

^"^ = E E Un) {0^{n)) , (1) 

c=l,8 n 

each factorized into two parts, fc{n) and {0^{n)), respectively, describing the production 
of a QQ pair in a state \nc) and the transition of this state into the quarkonium H. Here 
c = 1 (c = 8) denotes a colour- singlet (colour-octet) QQ pair. Besides spin S, orbital 
angular momentum L, and total spin J of the QQ pair, one more index, the dimension 
of the operator O, is required in order to specify n. The short- distance part fc{n) is 
calculable in perturbation theory in as{fi) with /i of the order of the heavy-quark mass 
m. For a non-relativistic system such as the heavy quarkonium, the heavy-quark velocity 
V (in the quarkonium rest frame) is the natural scale, which governs the magnitude of the 
long-distance matrix elements (MEs) {0^{n)). For calculations accurate to a given order 
in V, only a finite number of terms contribute to (|l]). 

The above-outlined factorization into the short-distance physics of the QQ pair and 
the long-distance physics of bound-state formation is a rigorous result of non-relativistic 
QCD (NRQCD) for inclusive quarkonium decays in the limit of large mass m of the heavy 
quark where the multipole expansion is valid[|T|. For the case of (inclusive) quarkonium 
production one needs the additional assumption |l], Q that all singularities, which neither 
cancel between real and virtual contributions nor can be absorbed into the MEs, can be 
factorized into the parton distribution functions (PDF). A study on an extension of the 
factorization approach to exclusive decays 0] has recently started. 

Two facts determine which MEs actually have to be included for an 0{v^) calculation 
of a given process. First the scaling of the MEs with v, and second the magnitude of the 
respective short- distance coefficients. The latter is determined by the powers of as{m) 
and, possibly, by kinematical factors. 

2 Short-distance processes 

Let us consider a few examples of short-distance reactions. The lowest-order short- 
distance processes relevant for quarkonium hadroproduction are 
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qq ^ QQi(X) X = 'S^,'D^,... 

QQ8(X) X = ■'S^,'D,,... al 

gg ^ QQi(X) X = 'So,'Po,2,'D2,... al 

QQ,iX) X = ^So,'Pu'Po,2,''^'Dj,... a2 



(2) 



The selection rules in (^) are easily understood. Annihilation of quarks and antiquarks 
proceeds through one-gluon (one-photon) exchange and is, hence, restricted to J"^ = 1^ 
states. Two gluons in a colour- singlet state are charge- conjugation (C)-even, while two 
gluons in a colour-octet state can be C-even or C-odd, depending on whether the d or f 
SU{3) structure constants are involved. Two-gluon production of the ^Pi state is forbid- 
den by virtue of the Landau- Yan theorem, independent of colour (provided the incident 
gluons can be taken on the mass shell!). The vanishing of the squared matrix element 
for gg -^ ^Si + g (again, true for on-mass-shell gluons) seems to be an accidental can- 
cellation between the two QED-like diagrams and the diagram involving the three-gluon 
coupling, which appears only in 1~ production. It would be interesting to check whether 
the rates for ^Di production and the O(w^) correction to ^5*1 are non-zero. Quarkonium 
production through ^Pi and D-wave intermediate states has not yet been included in 
phenomenological analyses. 

In photoproduction we encounter the following reactions up to and including the order 
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4^ 



7g ^ QQ8(X) 'So,'Po,2,'D, 

QQi(X)+g ■'S,,^P,,^Dj,... ela^^al 

QQ8(X)+g 25+iLj... e2^«em«2 

QQ8(X)+7 35i,lPi,3D^... (e2j«em)'«s 

7q -* QQ8(X) + q iS'o,=^Pj,iZ^2... e\a,^al 

QQ8(X) + q 3^1, 3^1... 



(3) 



3Q. 3n. .;2^ 



The selection rules for the reactions involving only one gluon (first and fourth ones) are 
the same as in QED. The second reaction is again restricted by the definite C property 
of the two-gluon colour- singlet state. Photon-quark-initiated reactions differ whether the 
photon couples to the heavy-quark line (oc Cq, where eg is the electric charge of the heavy 
quark in units of e) or the light-quark line (oc Cq) . 

In electron-positron annihilation via the exchange of an s-channel photon or Z'^-boson, 
the following final states can be reached: 
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ig^A? 
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. 'S,,'P^,'Dj,.. 


. - 


as^ 
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35i,iPi,3D^,.. 
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. - 


c^U 


QQ8(X) + qq 
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's. 


a^ln^e 



(4) 



The selection rules for one-photon exchange processes are the same as for the analogous 
photoproduction processes. The same holds true for the vector-current part of T? exchange 
(oc gy). The axial- vector current contribution has opposite C-transformation property to 
the vector- current contribution. Therefore any QQ(X) state can be reached through Z° 
exchange. I have also indicated the behaviour with ^ = (2m/ -^i)^, where a/s denotes the 
e"'"e~ cm. energy. 

One immediate consequence of @ is the value of A in the angular distribution, 
do"/dcos6' oc 1 + A cos^ 6', of (direct) J/-?/^ production with respect to the beam axis 
at large z = 2£^j/^/y^ and for mj/^ <^ ^/s <ti mzo (CLEO energies, ^/s ~ lOGeV). If 
colour-octet production is allowed (as is the case in both the colour-evaporation model 
(GEM) and in NRQCD, see below) then ^^o + g production dominates, resulting[|, Q in 
A = l. Explicitly!!], for r < 1 

ckr ^ 16A^ Uoi/XSo)} + 4 (Oren))} (l + cos=«) . (5) 

dcosp 9ms^ I m^ J ^ ^ 

On the other hand, if colour-singlet production dominates then (QQ)i('^Si)gg is the 
leading process, resulting inp A = —(1 — ^)/(l + ,^). This is easily understood: the major 
configuration in the colour- singlet model (GSM) has the two gluons recoiling against the 
J/ip with about equal energy. For ^ — > 0, i.e. neglecting the J/ip mass, conservation of 
parity and angular momentum then predicts A = —1. Hence, the measurement of the 
angular distribution appears as a gold-plated test of the colour-octet mechanism. 

However, also in the GSM differential cross section there is an A = +1 component, 
namely when the energy of one gluon is at its lower limit (and, necessarily, the other 
gluon's energy at its maximum value). Integrating the differential cross section over the 
region Ei < mv'^, where Ei are the gluon energies in the QQ rest frame, I find in the limit 
r <^ 1 and v^ <S 1 the NRQGD result (^ with the curly bracket replaced by 

81 TV 

Thus, if physics emphasizes the soft-gluon region the GSM result reproduces the NRQGD 
one with ^^(mf^) ~ 1. 

Alternatively, we can assume that gluons up to energies k ~ mv are important. Inte- 
grating the squared matrix element over the region Ei < mv, I again obtain a NRQGD- 
type result, but this time with the curly bracket replaced by a term proportional to 
as{rnv)v'^\R{0)\'^. Identifying as{mv) ~ v we now predict (Cg (^S'o)) to scale like v^ 
rather than like w^ w.r.t. {Oi {^Si)). 

3 Velocity-scaling rules 

The velocity scaling of the MEs {0^{n)) = (C»f (d,X)) in NRQGD is determined by the 
number of derivatives in the respective operators and the number of electric or magnetic 



dipole transitions between the QQ pair produced at short distances and the QQ pair in the 
asymptotic quarkonium. The velocity scahng for the leading MEs of the most prominent 
quarkonium states are listed in Table |I]. These leading MEs up to and including the 
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Table 1: Relative scaling of (Cf (X)). 



order v'^ are the ones used in current phenomenological estimates of production rates, 
in particular those MEs, which are enhanced by the short-distance coefficients. Note, 
however, that there are non-leading MEs that are less suppressed than some of the MEs 
given in Table |[ For example, the ME (Cf (rf = 8, ^5*1)) is suppressed by only f^ relative 
to the leading ME (Of (rf = 6,^S'i)) given in Table ^ it therefore gives rise to the first 
0{v) correction to 1 production. Even though it has the same short- distance factor as 
the leading colour-singlet ME, it may be favoured kinematically and, therefore, result in 
a sizeable contribution 0. 

In the limit of small velocity, the NRQCD MEs for S'-wave states and, hence, pre- 
dictions for their production (and decays) reduce to those of the CSM. In the CSM all 
colour-octet MEs are zero, while the colour-singlet ones are related to the quarkonium 
wave function at the origin 






<ip{nS) , 



{Or^-^\lO,'D,)) = 5^-^ \KM\' ■ (7) 



(Or^^"^\8, ^P,)) = (2 J + 1) -ii Kp 

ZTT 



The scaling of the wave function and its derivatives at r = with v, e.g. |i?5(0)p oc (mv)^, 
|-Rp(0)p oc (mv)^, hence determines the scahng of the MEs in the CSM as shown in 
Table |I|. _ 

Opposite to the scaling rules of the CSM are those of the CEM. Here any QQp(X) 
state has probability f " to reach any quarkonium state H. The only suppression that 
occurs is the one in the production of the QQ(,(^^+^Lj) state oc v^^, see Table |l|. 

The velocity scaling rules (VSR) of NRQCD were derived [Q on the basis of consistency 
requirements of the NRQCD Lagrangian. These are certainly valid for Aqcd "^ f^^v^ "^ 
mv, i.e. when as{mv) is small enough to renormalize the power ultraviolet-divergent 
contributions perturbatively. Alternative velocity scaling rules are possible, depending 
on the sizes of A and v. Here A is "the typical" scale of QCD, i.e. the one that sets 
the scale for higher-twist corrections, A may range from Aqcd to a typical hadronic 
scale and, hence, be as large as 1 GeV. I consider two cases. First, mv"^ falls below A, 
but A is still small with respect to mv, i.e. as{mv) <^ 1. Then the ordinary multipole 
expansion is still valid, but it is A that cuts off soft gluons rather than the binding 
energy^, leading to a double expansion in f ~ {p)/m and A = k/{p), where /c ~ A is the 
typical momentum of the dynamical (non-Coulombic) gluons. The scaling rules for this 
alternative scenario are also displayed in Table |l| (labelled VSRl). These can be argued 
in two ways. Either by considering the energy shift of the QQg Fock state a la BBL|[1| or 
from dimensional arguments. These tell us that electric and magnetic dipole transitions 
scale as r(ii^l) oc g^k^r^ and r(Ml) oc g'^k^/m'^, respectively. Hence the probabilities T/k 
scale as A^ for El and as A^f^ for Ml transitions owing to r~^ ~ mv. For k ~ m^v"^, the 
standard NRQCD scaling rules are reproduced. 

A second set of velocity-scaling rules (VSR2 in Table |1|) is arrived at when as{m,v) 
is not small enough to allow for a perturbative treatment. In non-relativistic QED, the 
power UV-divergent terms appearing in diagrams containing "soft" photons with energy 
k ~ mv can be cancelled perturbatively by corresponding terms in the bare coefficients^. 
This leaves only "ultrasoft" photons of energy k ~ mv"^ for which the usual counting rules 
of the multipole expansion hold, ~ v"^ for El and ~ f ^ for Ml transitions. However, for 
the counting rules of soft photons only factors of e and 1/m enter, so that ip'^{g/m)p- Aip 
and ip^g / m)a ■ Bif) contribute to the same order. Fock states that can be reached by both 
El and Ml transitions are suppressed by v^. Since nothing can be said about the velocity 
scaling of the non-perturbative renormalization required in NRQCD for as{mv) not small, 
y3 j-gjxiains the scaling of spin-flip transitions. Since the gA-piexjiv is not renormalized it 
still scales as f^.0 These general findings are in agreement with the explicit estimate of 

^ This possibility was suggested by Martin Bcnckc. 

^ After completion of this work I learned about an erratum to Ref.M, in which spin- flip transitions 



the scaling of (Og (^Sq)) at the end of section 2, where it was found that this ME scales 
as f '^ (v^) if gluons of energies mv"^ {mv) are important. 

One should bear in mind that the velocity is not a fixed quantity for a given quarko- 
nium system. The velocities of 1 and J"*""*" states are, in general, different and such is 
the case for different radial excitations. Universal velocity-scaling rules are only meaning- 
ful as long as velocity differences are small compared to the actual velocities. It may even 
happen that the relation between A and mv or mv"^ differs from one state to the other. 

In order to study these problems I calculate the velocity v, the momentum p and the 
kinetic energy T of the heavy quark for the lowest-lying S-, P-, and D-wave quarkonium 
states, for a class of potentials of the form V{r) = Xr'^ + Vq. A power-like potential 
successfully describes quarkonium spectroscopy with[0 u ~ 0.1. Uncertainties of the 
quark-potential model description can be investigated by studying the power-dependence 
within a reasonable range, say Au = ±0.2. For every value oi u I use the IS* leptonic 
width and the IS" and 25" masses to fix the three remaining parameters m, A, and Vq. Note 
that I normalize the leptonic width by the quarkonium mass rather than the heavy-quark 
mass and do not include the known 0{as) correction. 

The results given in Figs. ^ and ^ show that the velocities are rather uncertain: varying 
z/ between —0.1 and 0.3 changes, e.g. v"^,^ from 0.36 to 0.23. An additional uncertainty is 
caused by the 0{as) correction K = 1 — {16/3)as/n to the leptonic width, which enters 
as K'^^^ ~ 0.63 for charm and ~ 0.76 for bottom. It can also be seen that the difference in 
the 15* and 25* squared velocities is non- negligible, it can be as large as 20%. Moreover, for 
charmonia we observe f ^ ~ v/2 ~ asim). The kinetic energy T is about 350 MeV for both 
the charm and the bottom systems. For bottomonia, however, f ^ is small enough for the 
NRQCD relation v ~ asimv) to be fulfilled. Hence, bottomonia seem to be heavy enough 
for the NRQCD velocity-scaling rules to hold, while the charm-quark mass is possibly too 
light for the notion of a single velocity to make sense and a universal scaling to hold for 
all charmonium states. 

4 NRQCD matrix elements 

The long-distance MEs can either be calculated using non-perturbative methods or ex- 
tracted phenomenologically from data. In order to estimate the total 1 production 
rates we need also total production MEs, which include the feed down from higher states. 
For T(IS') the total MEs are defined as follows: 

^'^^^''^^'^Mtot = ^ {{Ol^'''\'Si))BT[T{nS) -. T(1S)X] 

ra=l 



have also been assigned the probabilty v^. I am indebted to M. Beneke for bringing this work to my 
attention. 



lT(lS),3 



{Ol^'^'i'S^)) 



tot 



^^"("^)e5i))Br[xbj(nP)-T(lS)X] 
J2 {Oj^'"^\^Si)) Br[T(nS) ^ T(1S) X] . 



j=o 



n=l 



The inclusive branching ratios needed for (^ and the other T{nS) states are compiled in 
Table 0. 



Xb2(3P) 
Xbi(3P) 

Xbo(3P) 

Xb2(2P) Xbi(2P) Xbo(2P) T(25) T(15) 

T(35) 11.4 ±0.8 11.3 ±0.6 5.4 ± 0.6 10.6 ± 0.8 11.7 ±0.5 

Xb2(2P) 16.2 ±2.4 12. ±1.3 

Xbi(2P) 21. ±4. 15. ±1.9 

Xbo(2P) 4.6 ±2.1 2.3 ±0.9 



T{3S) 


T{2S) 


T{1S) 


16.2 


5.4 


5.4 


21 


7 


7 


4.6 


1.5 


1.5 





Xb2(lP) 


Xbi(lP) 


Xbo(lP) 


T{IS) 


T{2S) 


6.6 ±0.9 


6.7 ±0.9 


4.3 ±1.0 


31.1 ±1.6 


Xb2(lP) 








22. ±4. 


Xbi(lP) 








35. ±8. 


Xbo(lP) 








<6. 



Table 2: Branching ratios (in per cent) in the bb system; sequential decays have been 



included for the T(3S'), T(2S'), and Xbj(2P) branching ratios; based on the PDG data[|ri 
and, for the Xbj(3P), on Ref.|T^. 



4.1 Non-perturbative methods 

Since production MEs cannot be determined in lattice gauge theory one has to resort 
to less rigorous non-perturbative models. The leading-order (in v) colour-singlet MEs 
can be estimated in the quark-potential model since they are simply related to the wave 
functions at the origin, see (|^). Tables |^ and ^ list the values obtained|T^ from a QCD- 
motivated potential [1^. Uncertainties in the wave functions arise from several sources: 



(0f/'^(35i)) 1.16 GeV^ (Cf/'^(35i))[^^ 1.59 GeV^ 
(Cf^(3Pi)) 0.32 GeV^ {Of^^\''Si)) 0.76 GeV^ 

Table 3: Colour- singlet long-distance matrix elements for charmonium production 



tot 



tot 



tot 



11.1 GeV^ 



{0'^^^^\^Si)) 9.28 GeV3 {0'^^^^\^Si)) 

^0Xbi(iP)(3p^)^ 6.09 GeV^ 

{0^''^^\^Si)) 4.63 GeV3 {Oj^^^\^Si)) 

^0Xbi(2P)(3p^^^ 7.10 GeV^ 

{0'^^^^\^Si)) 3.54 GeV3 {0^^''^\'Si)) 

(Cf"^^^^(3Pi)) 7.71 GeV^ 

Table 4: Colour-singlet long-distance matrix elements for bottomonium production 



5.01 GeV^ 



3.54 GeV^ 



relativistic corrections oc f^, the form of the potential, and badly determined input pa- 
rameters. Among the latter the 0{as) correction to the leptonic width gives by far the 
biggest uncertainty: the factor 1 — (16/3)a;s/7r can change the square of charmonium wave 
functions by as much as a factor of 2! 

In order to assess the size of the other effects I calculate the wave functions for a 
power-like potential V{r) oc Xr'^ and study them as functions of u. In this way we modify 
not only the form of the potential but also change the velocity v and, hence, the size 
of relativistic corrections. The IS* leptonic width and the 15* and 25* masses are used 
as input. The results for |i?„/(0)p = d'Rni(r)/dr^|r=o are shown in Figs. |l] and |^. In 
the case of charmonia, the nS to 15* ratio is rather stable but the ratio nP to IS* is 
quite uncertain, by a factor of almost 2. This uncertainty is not much smaller for the 
b system, in particular for the higher radial excitations. The D-wave wave functions 
are mostly uncertain. This underlines that the uncertainties of colour- singlet MEs must 
not be forgotten when estimating the ranges of theoretical predictions for quarkonium 
production. 

Colour-octet production MEs have so far been estimated only phenomenologically. 
Here I present the first non-perturbative estimates obtained in two different models. One 
calculation of the colour-octet ME, in collaboration with P. Hoodbhoy||15||, is based on 



the method of Leutwyler-Voloshin|TE 
gluon condensate (asG^). We obtain 



Here colour-octet MEs are given in terms of the 
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Table 5: Colour-octet long-distance matrix elements {O^i^^^^Lj)) for charmonium pro- 
duction in units of 10"^ GeV^+^'^. 



{Ol'^'Po)) 



12 



625 



m' {Ol'^CS,)) . 



(9) 



Note that (Og (^Pq)) is small and, as a consequence of (5^) = —{E'^), negative. 

Long-distance MEs can also be estimated by comparing NRQCD calculations with 
CEM ones. Leaving the details to a further publication I simply quote the results 
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(10) 



Here A = m{HQ) — mq is the difference between the mass of the heavy quark and that of 
the lightest meson containing it. The long-distance factors F[H] can be extracted from 
the ratios of the production rates for the various T{nS) {iIj{2S), J/^) states known from 
fixed-target experiments together with spin symmetry and an ansatz for direct T{nS) 
production, which relates these to the respective leptonic widths. The results are given 
in Tables ^ and ^ In the case of bottomonia, I present my estimates as a function of 
the size of the cross section of the as yet unobserved 3P Xbj states. The models CEMO, 
CEMl, and CEMIOO are defined by R^p = 0, 1, and 100, respectively, where 
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We consider CEMl as the most sensible choice. 
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Table 6: Colour-octet long-distance matrix elements {O^i^^^^Lj)) for bottomonium pro- 
duction in units of 10~^ GeV^"*"^^. 



4.2 Phenomenological determination of MEs 

Tables ^ and |^ quote also the values of colour-octet MEs determined in a study of quarko- 
nium production at the Tevatron|T^. The usual assumption rn? {O^ {} Sq)) = 3{Og(^Po)) 
is made and, for bottomonium, I have identified the 3L MEs with the corresponding 2L 
ones (L = S,P). 

The values of phenomenological determinations of (in particular colour-octet) MEs 
should be considered with great care. Numbers found in different processes may deviate by 
large factors and, yet, be consistent if proper error estimates were indicated. Uncertainties 
arise from several sources. First, a substantial uncertainty is caused by the large value of 
as{fJ,), which enters the various processes with different powers. Realistic error estimates 



should take into account the variations of Aqcd and the scale fx ~ Jm^ + p^. 
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A second source of uncertainty is introduced by the parton distribution functions 
(PDFs). Quark and gluon distribution functions enter in different combinations, at differ- 
ent scales, and in different x regions. Hence, comparisons of B decays into cfiarmonium 
(independent of PDFs), pliotoproduction of J/ip at low energies (dominated by single 
parton-distributions at large x), and J/ip production at the Tevatron (governed by the 
product of two gluon-distributions at small x) are highly non-trivial. 

Third, essentially all analyses are leading order in both v and as(m) but, the size 
of higher-order corrections in v and/or a^ may well be process-dependent. Moreover, it 
should be stressed that existing calculations are separately of leading order in v for the 
colour-singlet and colour-octet contributions, i.e. do not consistently include all terms of 
a given order in v. 

Last but not least, there exist higher-twist corrections, which enter at different levels. 
Examples are diffractive contributions to the z — > 1 behaviour of J/ip photoproduction 
oc A^/(l — z)ml, A/(l — XF)fni effects at large-x^ charmonium production, and effects of 
intrinsic kx at low-p^ quarkonium production at the Tevatron. 

5 Discussion 

We start by comparing the values of the charmonium colour-octet MEs extracted from 
the Tevatron with the non-perturbative results. The MEs {O^i^Si)) come out quite 
similar for H = J/ip, "^(25"), and Xcj- This is not surprising since ratios of production 
cross sections for these states H are rather universal|T^ and, hence, the MEs are well 



constrained. On the other hand, where there are less data, differences show up: the 
Tevatron values for {Og{^So, ^Pq)) are considerably larger than the non-perturbative ones. 

Turning now to bottomonia, we find similar values for {Oj{^So)) if R^p (see (0)) is 
close to our preferred value ~ 1. This also holds for the actual combination extracted 
from the Tevatron, (Cj(^S'o))/3 + {Oj{^Po))/ml, which was foundQ to be 7.9, to be 
compared with 16, 10, 0.31 (xlO~^GeV^) for R^p = 0, 1, 100, respectively. 

For the choice R^p = 1, however, the Tevatron value for {Og{^Si)) is much larger than 
the CEMl one. On the other hand, the situation is reversed for (Cj('^S'i))|dir so that, in 
fact, the values for (Cj(^S'i))|tot are rather similar. Again not a surprise, given that so 
far only total (direct plus indirect) T{nS) cross sections have been measured. Since the 
Tevatron values correspond to very large production cross sections for the as yet unob- 
served Xbj(3P) states, it is likely that the ^^i octet MEs for Xhj{nP) are overestimated 
while the ones for T{nS) are underestimated. 

Finally, I consider prompt quarkonium production from Z*^ decays at LEP and calculate 
branching ratios in both NRQCD and the GEM. I improve previous calculations in the 
CEM|g, |2g, |, ^ and the CSM[gg, ^ 0, § /NRQCD |2§ by carefully including in the 



3/ip and J2 ^ cross sections the feed down from higher states. Moreover, branching ratios 
for various other states are also presented. Also the GEM prediction for the (QQ)i + QQ 
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ZO^ 


(QQ)sS 


Z"^ 


(QQ)3qq 


ZO^Q 


(qq)^q 


Data 




CEM 


NRQCD 


CEM 


NRQCD 


CEM 


NRQCD 




J/M,^ 


0.014 


0.023 


13.6 


8.4 


0.14 


5.4 




Jli^L, 


0.023 


0.028 


22.7 


18 


0.23 


7.8 


19±10 H 


^(25) 


0.0033 


0.006 


3.3 


5.9 


0.033 


3.6 




XcO 


0.0046 


0.002 


4.5 


4.2 


0.046 


0.99 




Xcl 


0.014 


0.006 


13.6 


12.5 


0.14 


1.13 




Xc2 


0.023 


0.010 


22.7 


20.9 


0.23 


0.31 




T(15)|di, 


0.030 
0.020 
0.00059 


0.0076 


0.60 
0.39 
0.012 


0.034 


0.014 

0.0090 

0.00027 


0.76 




ET 


0.081 


0.14 


1.6 


4.8 


0.037 


1.7 


10 ± 5 [0 



Table 7: Branching ratios (in units of 10 ^) of 7P decays 



final state is calculated here for the first time. The results (Table |^) are obtained with the 
help of the MEs given in Tables l^-^. For direct T(IS') production, estimates corresponding 
to i?3P = 0, 1, and 100 are given. 

It can be seen that NRQCD and CEM predictions are rather similar for the colour- 
octet processes, but differ substantially for the colour- singlet one. Both predictions are 
compatible with the data. Concerning the NRQCD predictions: one can see that the ratio 
of colour-singlet to colour-octet production decreases when going from J /ip to the sum of 
the T states. This is in contrast with the expectation, since the octet contribution scales 
as v^ \t? ^ (\/^ = 2mQ/ ^/s) and both factors are larger for charmonium. This supports 
the observation made above that the colour-octet MEs of the P-wave bottomonium states 



as obtained [jT^ from the Tevatron are overestimated. 
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Figure 1: Average squared velocity, momentum, and kinetic energy of the charm quark 
as functions of v for charmonium states: S*- waves (sohd hues), P- waves (dashed hues), 
D-waves (dotted hues). Also shown is |i?y(0)p = d'Rni(r)/dr'|r=o- 
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Figure 2: Same as Fig. |l] but for bottomonia. 
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